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Abstract 

The Energetic Neutral Atom (ENA) full-sky maps obtained with the Interstellar Boundary Explorer 
(IBEX) show an unexpected bright narrow band of increased intensity. This so-called ENA ribbon 
results from charge exchange of interstellar neutral atoms with protons in the outer heliosphere 
or beyond. Amongst other hypotheses it has been argued that this ribbon may be related to a 
neutral density enhancement, or H-wave, in the local interstellar medium. Here we quantitatively 
demonstrate, on the basis of an analytical model of the principal large-scale heliospheric structure, 
that this scenario for the ribbon formation leads to results that are fully consistent with the observed 
location of the ribbon in the full-sky maps at all energies detected with high-energy sensor IBEX-Hi. 


1. INTRODUCTION 

The Interstellar Boundary Explorer (IBEX) has pro¬ 
vided the first energy-resolved all sky maps of the flux 
of energetic neutral atoms (ENAs). These IBEX maps 
reveal, above a general solar wind-structured ENA flux 
background, a ‘ribbon’ of increased flux. Several models 
have been proposed to explain the source, location, and 
structure of the ribbon. 


McComas et al. (2009a) and Schwadron et al. (2009) 
were the first to suggest that the ribbon might result frorn 
consecutive charge-exchange processes. This scenario 
has, subsequently, been modelled quantitatively (see. 


e.g., Heerikhuisen et al.| 

2010| IHeerikhuisen & Pogorelov 

2011 IStrumi 

k et al.|20i; 

[[IMobius et al.|20131|Schwadron 


& McComas 


.120131 IZirnstein et al.||20131 IBurlaga et ah, 

2014| Heerikhuisen et al. 2014| Isenberg| 2014). Despite 
the modehs basic success to explain the EN A ribbon, 
there remain critical open questions regarding the sta¬ 
bility of the (pick-up ion) seed distributi on of the ribbon 


ENAs in the local interstellar medium (Elorinski et al 


2010 Gamayunov et al. 2010 Burlaga JNess 


2(JTT) 


Other suggestions comprise the ideas that the source re¬ 
gions of the ribbon ENAs are located far beyond the he- 
l iopanse (HP) at the edg e of the local interstellar cloud 


(Grzedzielski et al. 


(kahr et al. 


mry 


20101) or rather inside the heliosphere 


Kucharek et al. 2013 Siewert et al. 


z013[ ). A diherent scenario that involves the magnetic 


and neutral density structure of the local interstellar 
medium (LISM) on the one hand, but assumes the pro¬ 
duction region of the ENAs to be mainly in the inner 
helioshea th (IHS) on the other hand was suggested re¬ 
cently by Eichtner et al. (2014). 

All of these hypotheses have been summarize d and crit¬ 
ically assessed in d etail in the review papers by McCom^ 
et al.| (|2014a|) and jMcComas et al.| ( 2014b|) with the re¬ 
sult that, while it is clear that the ENAs establishing 
the IBEX ribbon are related to the interaction of the he¬ 
liosphere with the local interstellar medium, there is no 
consensus yet on their source region(s). 

In this paper we follow up on the idea that the IBEX 
ENA ribbon is a result of a so-called H-wave (section 


2) tr ansiting through the heliosphere (Eichtner et al 


2014). In order to translate that idea into a quantitative 


model, we construct the geometry of the ribbon (sec¬ 
tion 3) within the framework of a simple but well-suited 
model of the principal large-scale heliospheric structure. 
We discuss a best fit to the ENA data obtained with 
the IBEX-Hi detector along with the significance for the 
findings regarding the relation between the ribbon and 
the local interstellar magnetic field (section 4) and the 
sensitivity of the results to parameter changes (section 
5). Section 6 contains a brief summary of results and the 
conclusions regarding the H-wave hypothesis. 

2. MODEL OF THE HELIOSHEATH AND THE H-WAVE 
INDUCED RIBBON FORMATION 

2.1. Plasma flow 

For the present purpose both the solar wind plasma 
flow in the IHS and the interstellar plasma flow in the 
outer heliosheath (OHS), i.e. in the regions between the 
termination shock (TS) and the heliopause (HP) and out¬ 
side the latter (formally up to the interstellar bow shock), 
can be described as being incompressible (V • u = 0, see 



Figure 1. The flow lines in the inner and outer heliosheath (see 
Eq. 4). The outer thick black line is the heliopause (77 = 1, see 
Eq. 5). The red color corresponds to inner heliosheath (77 < 1) and 
the blue color to outer heliosheath (rj > 1). The inner circle is the 
termination shock. 
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Roken et al.||2Q15 ): 


u = -V* = -l|)e,- 


kpULISM 


dxs 
-(r^ - kx3)ey, 


( 1 ) 

( 2 ) 


with the scalar velocity potential in cylindrical coordi¬ 
nates 


^{p,X 3 ) = ulismXz + ■ k = const. (3) 

r 

where r = ^ p = and {xi^X 2 ^x^) 

denote Cartesian coordinates and ulism is the speed of 
the undisturbed LISM flow. In this formulation kuLisM 
is interpreted as the speed of the shocked solar wind in 
the IHS. The resulting flow lines x^{p) can be obtained 
from the equation 


dx2, _ 

dp Up 


or from the associated stream function (see Appendix A) 
with the solution: 



(4) 


The parameter p is identifying the flow lines and k char¬ 
acterizes the relative strength of the solar and the inter- 

i 


stella r wind, which we take as k = 2 (see [Roken et ah 


2015). With this and 77 = 1 we obtain the following for¬ 
mula for the HP surface 


1 



Figure 1 illustrates the resulting flow lines in the IHS and 
the OHS. The interstellar flow comes from the positive ^- 
direction. The black lines indicate the HP p = 1 (see Eq. 
5) and the TS, respectively. For simplicity, we assume 
the latter to be a Sun-centered sphere. The red lines in 
the IHS correspond to 77 < 1 , the blue ones in the OHS 
to 77 > 1 (see Eq. 4). 


2 . 2 . Structure in the neutral gas 


Eichtner et al. (2014) have summarized the arguments 
supporting the ^lew that the LISM is inhomogeneous 
and is likely to exhibit propagating wave- or pulse-like 
structures. These authors argued that the waves in the 
plasma must be expected to induce (via charge exchange 
cou pling) associated struct ures in the neutral gas (see 
also Shaikh & Zank 2010). In particular, as a conse- 
quence of a slow wave m the plasma (that propagates 
along the magnetic field tha t is oriented as sketched in 
the top panel of Eig. 2 , see McComas et al.j (2009b)), 
there should be a wave in the neutral gas, too, a so-called 
H -wave. The bottom panel of Eig. 2 shows a sketch of 
such an H-wave, representing an enhancement of inter¬ 
stellar gas density that propagates along the interstel¬ 
lar magnetic field and, after decoupling from the plasma 
that flows around the ‘obstacle’ heliosphere, penetrates 
the heliosphere. While in this sketch the H-wave front 
is perpendicular to the undisturbed interstellar magnetic 
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Figure 2. Top panel: The flux of energetic neutral atoms pro¬ 
jected onto the heliopause around which the interstellar magnetic 
held lines are draping (adopted from McComas et al. 2009b). Bot¬ 
tom panel: Sketch of the H-wave scenario in a plane perpendicular 
to the orientation of the undisturbed interstellar magnetic held 
(IMF, black lines). An enhancement of interstellar density, i.e. an 
H-wave of thickness a is propagating through the heliosphere that 
is depicted by the terminaton shock (TS) and the heliopause (HP). 
In the intersection region, which is indicated by the shaded areas 
and which forms a ring-like structure in 3D, the production rate of 
ENAs is increased due to the higher neutral density. 


field, not only this but also other orientations are quan¬ 
titatively studied below. Th e production rate of ENAs is 
highest in the IHS (see, e.g. [Sternal et al.|| 2 Q Q^ and di¬ 
rectly proportional to the neutral density (e.g. Fahr et al. 
2007) so that one must expect an increased ElNA flux to 
be generated in the shaded regions shown in the bottom 
panel of Eig. 2. 


3. DESCRIPTION OF THE RIBBON GEOMETRY 
Eigure 3 illustrates the b asic geometry o f the r ibbon 


in the scenario suggested by Eichtner et al. 
TS is assumed as a Sun-centered sphere wit 


([Mil). 

bn radii 


The 
ms R, 


the HP as an axisymmetric surface defined by Eq. (5), 
and the H-wave is indicated as the thick black line E. 
The angle between the directions to the intersection of 
E with the TS and the normal vector n is denoted by 
a, defining the ’inner boundary’ of the ribbon in a given 
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Figure 3. Geometry of the intersection of the H-wave with the 
heliosphere: In the chosen plane the termination shock (TS) is the 
inner circle, the heliopause (HP) is the outer curve and the H-wave 
is indicated by the plane E with a normal vector n. The angles 
a, 7 , and 7 ' represent respectively the ’inner boundary’ and ’outer 
boundary’ of the ribbon and, thus, its angular width. For more 
details see text. 


plane: 


a = arccos 



(6) 


Here d denotes the shortest heliocentric distance to the 
plane E. The ’outer boundary’ of the ribbon is defined 
by the directions to the intersections of the plane E with 
the HP. In a given plane those directions can be spec¬ 
ified in terms of the angles 7 and 7 '. If yhp and r'^p 
denote the heliocentric directions to these intersections 
the corresponding angles follow from 


7 ^') = arccos 


HP ■n 

|r(')^p||n| 


(7) 


Given the chosen symmetrical TS, the angle a is the 
same for both ’sides’ of the ribbon in a given plane, while 
the angles 7 and 7 ' are different. With a, 7 and 7 ' it is 
now straightforward to formulate the condition that a 
given line of sight intersects the ribbon: 

a</?<7>7' (8) 


where (3 is the angle between the chosen heliocentric di¬ 
rection r and the normal n (see Eig. 3), i.e. 


P = arccos 


Xifli + X2n2 + XsTls 

\J {xl + + xl){nl + n| + n|) 

where in the Gartesian coordinates as defined in sec¬ 
tion 2 . 1 : 

n= ( 77 . 1 , 77 . 2 ,^ 3 ) and r = {xi,X 2 ,x^) ( 10 ) 

With these formulas we are now in the position to cal¬ 
culate the location and angular width of the ribbon in 
the all-sky ENA flux maps as observed with the IBEX 
spacecraft. 



= arccos 


r • n 



Figure 4. The observed ENA ribbon in an all-sky map of the 
observed ENAs fluxes (ENAs/(cm^ s sr keV) at 0.71 keV vs. dif¬ 
ferent model results: Top panel: The ribbon (two thick black lines) 
as resulting from an inflnitely thin H-wave with (7 = 0, i.e. a plane 
E with d = O.IR (see Fig. 3). Note that also for this case the 
ribbon has a flnite angular width because the angles a, 7 , and 7 ' 
are different. Middle panel: The ribbon resulting from an H-wave 
with flnite thickness a = d 2 — di = OAR — O.IR = 0.3R is the re¬ 
gion between the black lines indicating its inner (upper) and outer 
(lower) boundary. Bottom Panel: The ribbon location and thick¬ 
ness resulting from a transformation of the heliopause deflned with 
Eq.(5): While the latter results in a lower boundary indicated by 
the green line, one obtains the red and blue lines when it is tilted 
by 5 and 10 degrees, respectively. The untilted but polynomial 
heliopause (see text) leads to almost identical result (black line) 
than Eq.(5), i.e. the green line. For these and all following all-sky 
maps, the data of which is available via the IBEX data release web¬ 
site http://ibex.swri.edu/researchers/publicdata.shtml an Aitoff- 
projection was used such that the heliospheric nose direction is in 
the center of the map (see Appendix B). 


4. THE RESULTING RIBBON GEOMETRY 

Eor the visualisation of the resulting band of higher 
ENA fluxes in the all-sky maps we use an Aitoff pro¬ 
jection (see Appendix B) for all following figures. As a 
first step we check on the principal location and width 
of the ribbon originating from an H-wave, i.e. we choose 
its thickness a = d 2 — di (with di ^2 denoting the shortest 
heliocentric distances of its sunward and anti-sunward 
side, respectively) and its orientation given by the nor¬ 
mal vector n. 

The top panel of the Eig. 4 gives the result of an in¬ 
finitely thin H-wave (i.e. with vanishing thickness a) ori- 
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ented perpendicularly to the vector n = ( — 1.8,1.3,1.5), 
which is corresponding to the unit vector n/|n| ^ 

(—0.672,0.485,0.560) and which is anti-parallel to the 
most likely direction of the undisturbed local interstellar 
magnetic field (see the discussion in section 5 below). Its 
shortest heliocentric distance is d = O.li?. Note that al¬ 
though the H-wave is infinitely thin, the corresponding 
ribbon is not: In an all-sky map it has a finite angular 
width because the angles a, 7 , and 7 ' are different, as 
is evident from the sketch in Fig. 3. The middle panel 
shows the result for an H-wave with the same orienta¬ 
tion but a finite width (as sketched in Fig. 2), namely 
a = d 2 — di = OAR — O.IR = 0.3R. The inner (upper in 
the all-sky map) angular ribbon boundary is determined 
by the angle a resulting from the intersection of the anti- 
sunward plane with the TS, and the outer (lower) bound¬ 
ary by the angles 7 , 7 ' resulting from the intersection of 
the sunward plane with the HP, as described in section 3. 
The bottom panel shows the effect of different transfor¬ 
mations of the HP on the lower ribbon boundary. First, 
a transformation to a polynomial shape is achieved by a 
Taylor expansion of the HP function defined with Eq.(5), 
for details see Appendix C. Comparing the green line (re¬ 
sulting from Eq.(5)) and the black line (Taylor-expanded 
HP function) reveals that the effect of a polynomial HP 
on the ribbon location and width is negligible. This can 
easily be understood as a consequence of the fact that, in 
the upwind heliosphere ( 2 : > 0 ), even in lowest (second) 
order the Taylor expansion is well approximating the HP 
defined with Eq.(5), see Eigure 7 in appendix C. Second, 
a tilt of the HP by 5 or 10 degrees, achieved by subse¬ 
quent rotations about the x- and ^-axis, results in the 
lower boundary indicated by the red and blue line, re¬ 
spectively. These_til^ngl^are motivated by the recent 
finding by Wood et al. (2014) that the heliotail direction 
is deviating at most by 10 degrees from the inflow direc¬ 
tion of the LISM. Evidently, the effect is significant and 
one concludes (i) from a comparison of the red and the 
green line that a tilt of the HP improves the agreement 
between the modelled and the observed ribbon and (ii) 
from a comparison of the red and the blue line that the 
tilt must be expected to be less than 10 de gree, consis¬ 
tent with the findings by Wood et al. (2014). Therefore, 
for all following computations, we used the HP function 
Eq.(5) with an additional tilt by 5 degrees. 

The results shown in Eig. 4 make it evident that the H- 
wave hypothesis results in an ENA ribbon at the correct 
location in the all-sky maps and that its width must be 


a few tens o f AU, as was already speculated in Eichtner 


et al. (2014). The given geometrical H-wave parameters 


were iterated such that the resulting ribbon geometry is 
a simultaneous best fit to the high-energy ENA maps at 
1.11, 1.74, and 2.73 keV provided by the IBEX-Hi sensor. 
The result for these are shown in Eig. 5, which demon¬ 
strates that the H-wave induced ENA ribbon geometry 
is generally consistent with the IBEX measurements at 
all these energies. Note that for the highest energy chan¬ 
nel at 4.3 keV shown in the bottom panel the ribbon 
feature is clearly present and well-fitted at high north¬ 
ern latitudes (top of the map), but that there are sig¬ 
nificant ENA emissions outside the ribbon at lower lati¬ 
tudes. These additional signals at higher energies are also 
known from m easurements with th e INCA instrument 
aboard Cassini ( Krimigis et al.| 2 Q^ ). An explanation of 



Figure 5. All-sky maps of the ENA fluxes (ENAs/(cm^ s sr keV) 
as observed by IBEX-Hi and the simulated best-fit (see text) ribbon 
geometry for the energies 1.11 keV, 1.74 keV, 2.73 keV, and 4.29 
keV from top to bottom. 


the associated ‘broadening’ of the ribbon at higher ener¬ 
gies probably requires invoking additional ENA sources, 
like secondl y ENAs as, e.g., discussed in [Heerikhuisen 


et al. (2014). 

TEe above findings corroborate the assumption that 
the ENAs forming the ribbon in the IBEX-Hi all-sky 
maps can indeed originate in the intersection region of 
an H-wave with the inner heliosheath. While this repre¬ 
sents a rather different explanation from all other scenar¬ 
ios that have been suggested, so far, it does confirm the 
relation of the ribbon to the local interstellar magnetic 
field that is needed in most other scenarios, too. One 
must distinguish, however, in the present model between 
the field direction an d the normal vector t o the H-wave 
front: As discussed in Eichtner et al. (2014) it is, in prin¬ 
ciple, possible that the H-wave front is not perpendicular 
to the magnetic field. The effect of a different orienta¬ 
tion along with different widths of the neutral density 
enhancement is discussed in the next section. 


5. SENSITIVITY TO THE H-WAVE PARAMETERS 

To illustrate the sensitivity of the result to the orien¬ 
tation of the H-wave front, i.e. to check on the goodness 
of the best fit, the upper two panels of Eig. 6 give the 
ribbon geometry for a wave normal vector n whose direc¬ 
tion differs by 5 and 10 degree, respectively, from the best 
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Figure 6. All-sky maps of observed ENA fluxes (ENAs/(cm^ s sr keV) at 1.1 keV and overlayed simulated ribbon geometry. The upper 
two plots show the result if the normal vector of the H-wave front is changed by 5 and 10 degrees from the best fit direction. The lower 
two plots display the results for the best fit direction but with the two different H-wave widths a = d 2 — di = 0.3R — O.IR = 0.2R and 
a = d2-di= 0A5R - O.IR = 0.35. 


fit direction. While the first result (upper left panel) is 
still compatible with the observational data, the second 
(upper right panel) is clearly not. 

It interesting to note that the best fit normal vec¬ 
tor and the upwind direction —ulism have an angle of 
(55 ± 5)^, which is close to the 49^ between the upwind 
direction for the untilt ed HP and the interstellar mag¬ 
netic field estimated by Heerikhuisen et al. (2014). The 
ecliptic longitude Xed ^ 205"" and latitude ded ~ 35*^ 


of the best fit normal vector are also slightly di fferent 
than the values discussed in the literature (e.g., Witte 



however, the simplifying assumptions regarding the TS 
and HP surfaces made in our analytical approach, one 
can safely state that the normal vector is closely related 
to the direction of the undisturbed local interstellar mag¬ 
netic field. 
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The two lower plots in Fig. 6 show the influence of 
the H-wave width and reveal that the region of increased 
neutral density as well as, in turn, of enhanced ENA 
production and flux must indeed be assumed to be about 
25-28 AU wide, assumin g a TS radius of 84 to 94 AU 
dStone et al.]2^poM ). 


6 . CONCLUSION 

In this paper we have constructed an analytic model 
that reproduces the correct geometry of the IBEX ribbon 
in the all-sky ENA flux maps and, thereby, corroborates 
the hypothesis that a propagating localized density in¬ 
crease in the neutral interstellar gas, termed an H-wave, 
can be the cause of the IBEX ribbon. The best fit of 
this geometry to IBEX ENA data depends particularly 
on the orientation and width of the H-wave whose transit 
through the heliosphere leads to an increased production 
of ENAs in the inner heliosheath. Despite the simplify¬ 
ing assumptions regarding the termination shock and he- 
liopause surfaces, the proposed scenario, although rather 
different from all others that have been invoked to ex¬ 
plain the IBEX ribbon, makes it likely that the ribbon is 
closely related to the direction of the undisturbed local 
interstellar magnetic field. 

In subsequent work we will extend the modelling to 
a computation of the actual fluxes of ENAs resulting 
from an H-wave intersecting the inner heliosheath, i.e. 
we will evaluate the relevant line-of-sight-in tegrals (see, 
Eichtner et al. 12014) from an in¬ 


e.g.,|Sternal et al. 120081 |Fi 
ner boundary (IB) at tne ( 


boundary (IB) at the detector to an outer boundary 
(OB) sufficiently beyond the ENA source region (i.e. in 
the present case beyond the HP): 

OB 

^{EeNA,'&,^) = ['rt-pfp{'>’p)'rt-H<yexVrel]ds ( 11 ) 

IB 

with the solar wind and pick-up ion proton velocity 
{vp) distribution function fp and number density n^, the 
charge exchange cross section aex and the relative speed 
Vrei between a proton and an interstellar neutral hydro¬ 
gen atom. There are two key ingredients: first, the pro¬ 
ton velocity distribution function whose evolut ion has 
to be computed fr om a transp o rt mo del like in Eahr fc 
Eichtner (2011) or Eahr et al. (2014), but here for the 
IHb. I'he structure of fp in the intersection region of the 
H -wave and the IHS determines the ‘fine structure’ of the 
ENA fluxes in the ribbon in the all-sky-maps at different 
energies. Second, the ENA flux is directly proportional 
to the number density uh of the hydrogen atoms and, 
thus, to the H-wave signature. Erom a s tronomical o bser 


vations (e.g. Haverkorn & Goss 20071 Welty |2Q07|) a nd 
corresponding simulations ( |Henne belie Auclit||2007 ) it 


is derived that uh can easily vary by a factor of two and 
more down to the few AU scale in the so-called warm neu¬ 
tral (interstellar) medium. Since the latter reflects the 
properties of the n eutral component in the LISM (e.g., 
Stanimirovic||2009|) , one can expect to see such variation 
as a local H-wave with a two- to threefold enhanced uh- 
This directly translates, via Eq.(ll), into a correspond¬ 
ingly increased ENA flux, which m turn represents the 
general ribbon feature in the all-sky maps. 
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APPENDIX 

A: ANALYTICAL CALCULATION OF THE FLOW LINES 
In cylindrical coordinates (p = X 2 ^xs) the velocity potential (Eq. 3) reads 

kULISM 


^{p,xs) = ulismXs + 




(Al) 


with ulism denoting the velocity of the undisturbed local interstellar medium (LISM) and k the relative strength of 
the (shocked) solar and interstellar wind. The flow lines are on surfaces of constant stream function T. The latter can 
be derived from the incompressibility condition V • u = 0 , which reads in cylindrical coordinates explicitly: 


Erom this it follows that the stream function must fulfil the two equations: 


which have the solutions: 


T = 


/ 


1 9T 

pdxs 

kULISMp^ 


= Uo 


1 5T 

P dp 


— Ha: 


<i*s + G(,) = -*5i|S^ + 0 {rt 


(A2) 

(A3) 

(A4) 


and 


T = 


= - [ ULisMpdp^ f -J^=^£^dp + F{x3) = --Ulism p‘^ - ^=^^+F(x 3 ) (A5) 

J J VP Fxi ^ x/p^ + xi 


^3 “ \/F ^3 

with two functions G{p) and F{x 3 ) occuring as integration constants regarding integration w.r.t. x^ and p, respectively. 
Chosing G{p) = \ulismP^ and F^x^) = 0 leads to 


. N kULISMXs 1 2 

T(p,X3) =- - -UlismP 

VP + ^3 ^ 

The condition T(p, xs) = f) = const describes the flow lines. With the definition p = —fj/{kuLisM) one has 

^2 




X3 


+: 

from which one finds the desired equation (4) for the flow lines 


(A6) 


(A7) 


X3{P)= )P( 1 - 


P 


2k 


(AS) 


B: COORDINATE TRANSFORMATION AND AITOFF PROJECTION 

The Cartesian coordinates introduced in section 2 are related via 

Ti = RsinP cos (p ; T 2 = R sin P sin p ; X 3 = RcosP (Bl) 

to a heliocentric spherical polar coordinate system. In order to plot the desired all-sky maps centered on the heliospheric 
nose the following transformation is applied 

x[ = X 3 = RsinP'cosifi' ; X 2 = X 2 = R sin P'simp' ; x'^ =—xi = RcosP' (B 2 ) 

with the new latitude and longitude angles P' and p'. 

The Aitoff projection maps these spherical polar coordinates on Cartesian ones (t, y) in a plane via 


2 (acos( 0 ) sin (|) 


y = 


M) 


a = arccos cos 0 cos — 


W' 


sin((a) ’ sin((a) 

The dependence of the angles (j) and A on P' and p' as well as on P and p is given by 

TT 

(j) = — — P' = arcsin(— cos p sin P) ; A = p' = arccos (cos P/ cos (j)) 
and they must be interpreted as the latitude and longitude from the central meridian, respectively. 


(B3) 

(B4) 



Sylla and Fichtner 


C: POLYNOMIAL HELIOPAUSE EUNCTION 
The Taylor expansion of the HP function Eq.(5) up to sixth order reads: 


Xz{p) = V2 - 


3^2 2 


5V2 4 


ly/2 

81 ^ 




(Cl) 


The corresponding ‘polynomial’ HP surfaces are illustrated with Fig. 7 where the parabolic HP (up to second order, 
green line) is plotted in the x — z-plane along with the expansion to fourth order (blue line) in comparison with the 
HP according to Eq.(5) shown as the red line. Evidently, even in lowest order the approximation is already reasonable 



0 

Figure 7. The HP curve defined with Eq.(5) (red line), its Taylor approximations up to second (green line) and fourth (blue line) order 
and an extremely narrow HP. 


in the upwind heliosphere, i.e. for positive 2 ). By multiplying one or more Taylor coefficients with factors greater than 
unity a narrower HP can be obtained, like the example illustrated with the violet line (obtained by multiplying the 
coefficient of the sixth order by 100). These HP shapes are not only strongly deviating from the HP function Eq.(5) 
and are, thus, strongly inconsistent with the flow field, but they still, like all polynomial surfaces, are characterized by 
a diverging cross section of the heliotail, an undesired feature that the HP defined with Eq.(5) avoids. 









